Weakly frustrated two-dimensional Heisenberg antiferromagnets: thermodynamic 

properties from a non-perturbative approach 
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We analyze the thermodynamic properties of the spin-S two-dimensional quantum Heisenberg an- 
tiferromagnet on a square lattice with nearest and next-nearest neighbor couplings in the Neel phase 
{J2/J1 < 0.4) employing the quantum hierarchical reference theory (QHRT), a non-perturbative 
implementation of the renormalization group method to quantum systems. We investigate the stag- 
gered susceptibility, the structure factors and the correlation length at finite temperature and for 
different values of the frustration ratio. From the finite temperature results, we also extrapolate 
ground state properties, such as spin stiffness and spontaneous staggered magnetization, providing 
an estimate of the extent of quantum corrections. The behavior of these quantities as a function 
of frustration may provide some hint on the breakdown of the Neel phase at zero temperature for 
larger values of J2. 

PACS numbers: 05.10.Cc75.10. Jm 05.70.Jk 75.40.Cx 



I. INTRODUCTION 

Several new magnetic materials characterized by com- 
peting (i.e. frustrating) interactions have been recently 
synthesizediiSi. Unfortunately, a detailed theoretical un- 
derstanding of the ground state properties of simple mod- 
els displaying both frustration and quantum fluctuations 
is still missing'^ and the interpretation of the experimen- 
tal results often relies on simple perturbative or semi- 
classical estimates. Needless to say, the effects of thermal 
fluctuations in these systems is far from being understood 
either, despite the fact that most of the experimentally 
accessible quantities explicitly refer to the temperature 
response. 

Among the different magnetic models with competing 
interactions which have been theoretically investigated, 
the so-called Ji — J2 model provides a paradigm of the 
open problems and challenging questions arising in this 
area of theoretical magnetism^. The Hamiltonian of the 
square lattice Heisenberg antiferromagnet with nearest- 
neighbors (-n.n) and next-nearest neighbors interactions 
(n.n.n.) is 



H, 



(1) 



where Ji, J2 > and S are spin-S" operators. The model 
is characterized by a single parameter: a = J2I Ji named 
frustration ratio. The interest in this model was trig- 
gered by the discovery of high-Tc superconductors and by 
the possible connections between the physics of the two- 
dimensional Hubbard model and the Cu oxides materials. 
It has been argued that the motion of holes in a quantum 
antiferromagnet gives rise to effective frustrating cou- 
plings in the undoped Heisenberg Hamiltonian"*, which 
eventually leads to the breakdown of Neel order. More 
recently, the physics of the Ji — J2 model raised a renewed 



interest due to the discovery of three vanadate com- 
pounds (Li2VOSi04, Li2VOGe04 and VOM0O4) which 
appear to be characterized by antiferromagnetic nearest 
and next-nearest neighbor interactions between spin-1/2 
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At the classical level, the zero temperature phase di- 
agram of the Ji — J2 model separates into two regions: 
for a < 1/2 the system shows Neel order with magnetic 
wave vector g — {tt,tt), while for a > 1/2 the two sub- 
lattice decouple, each displaying Neel order, leading to 
(infinitely) degenerate classical ground states. This de- 
generacy is lifted by quantum fluctuations which, at the 
semi-classical level, select collinear order with periodic- 
ity characterized by wave vectors (tt, 0) or (0,7r). In re- 
cent years the search for a new disordered phase in the 
intermediate region of frustration (a ~ 0.5), has been 
the main subject of theoretical investigation^ while the 
thermal properties of frustrated quantum models have 
attracted comparatively less attention. Finite tempera- 
ture studies cannot rely on the Monte Carlo simulations, 
successfully employed in the non-frustrated case: in fact 
the use of these methods is strongly limited by the well- 
known sign problem, a numerical instability which oc- 
curs both in fcrmionic systems and in frustrated bosonic 
models,-. On the other hand, the Lanczos-based ap- 
proaches do not suffer from the sign problem, but the 
limitations in the cluster size make rather difficult the 
calculations at finite temperature^. Nevertheless ther- 
modynamic parameters such as spin stiffness and spin 
velocity have been estimated by exact diagonalization on 
finite clusters followed by extrapolation to the thermo- 
dynamic limit. Unfortunately these small cluster results 
are always affected by the uncertainties in the finite- 
size scaling procedureSiifliii. Among the numerical ap- 
proaches, the high temperature expansion method seems 
to be the natural candidate for the analysis of thermal 
fiuctuations^. Indeed, the physics of the collinear phase 
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has been recently investigated by several high tempera- 
ture expansion studiesi2ii4. Besides the interest related 
to the physics of the vanadate materialsSifiii^, the study 
of the large a-limit of the model was also motivated by 
the possible observation of a finite temperature 'Ising' 
transition in the Ji — J2 model. In fact the residual 
Z2 symmetry of the Hamiltonian in the coUinear phase 
(due to the twofold degeneracy of the ground state) can 
be broken at finite temperaturei^, without violating the 
Mermin- Wagner theorem. 

On the analytical side, Ivanov and Ivanov^® inves- 
tigated the low-temperature thermodynamics for < 
a < 1, by a modified spin- wave theory (SWT). Us- 
ing Schwinger-boson mean-field theory Mila et al. ex- 
amined the spin dynamics of the model, as an effec- 
tive theory for the high-Tp Cu oxidesii. A comprehen- 
sive study of the thermodynamics of the collincar phase, 
within a (semiclassical) effective Hamiltonian approach 
(PQSCHA), has been performed by Capriotti et ai, also 
for the extreme quantum case S — 1/3^^, while only a few 
investigations have been devoted to the thermodynam- 
ics of the weakly frustrated model, in the Neel phase. 
In the low temperature regime some hint comes from 
the renormalization group (RG) formalism applied to the 
non-linear sigma model (NLSM). Within the NLSM ap- 
proach, the frustrating interaction is treated in mean field 
approximation and the fluctuations induced by the frus- 
tration are in fact neglected. The terms in the Hamil- 
tonian acting on the variables of the same sub-lattice, 
such as the J2 interaction, are integrated out in order 
to obtain the effective action and then the effect of the 
frustration reduces to a renormalization of the coupling 
constant g^. As a result, the inclusion of a next near- 
est neighbor interaction J2 > increases the bare cou- 
pling g and , following the RG picture, enhances short 
range fluctuations. Indeed, if the frustration is not strong 
enough to destroy the long range order, the low temper- 
ature physics of the systems is believed to be the same of 
the non-frustrated case, with physical properties, such as 
spin stiffness and spin velocity, renormalized by the inclu- 
sion of frustration. Spin- waves expansion results seems to 
support this picture. However, even in the non frustrated 
case, the picture based on NLSM and RG approach does 
not adequately represent the experimental behavior for 
the quasi-2D antifferomagnet with S > 1/2. Actually the 
agreement between theory and experiment seems to hold 
only for S = 1/2 in the low and intermediate temperature 
regions. In fact, the effective Hamiltonian description is 
generally valid in a regime characterized by strong short 
range correlations and it is not suitable to give a detailed 
account of the physics in a range of temperature where 
local order weakens or even disappears. 

In order to describe the interplay between frustration, 
temperature and quantum fluctuations, we follow here 
a different method which aims to reconcile the RG ap- 
proach with a microscopic description of the model. This 
theory, named Quantum Hierarchical Reference Theory 
(QHRT)^r does not rely on a coarse graining procedure 



leading to some effective action but instead is a way to 
implement the momentum space renormalization method 
directly to the physical Hamiltonian. Starting from the 
mean field approximation, we recursively include the fluc- 
tuations over all length scales by means of an evolution 
equation, which describes how the free energy of the sys- 
tem changes when fluctuations over growing wavelengths 
are taken into account. The original method developed 
by Wilson has been largely employed for the investiga- 
tion of the properties of quantum systems, but, contrary 
to the weak coupling RG method applied to effective 
field theories, our approach does not require any mapping 
on coarse grained actions and applies to all temperature 
regimes. The resulting evolution equation is formally ex- 
act but it is not written in a closed form: An ansatz 
for the momentum and frequency dependence of the dy- 
namical structure factor is needed in order to put the 
equation in a closed and solvable form. The choice of a 
parameterization for the spin-spin correlation function is 
the only approximation present in this approach. 

The paper is organized as follows. In Section // we 
derive the evolution equation for the Ji — J2 model. In 
Sec. /// we investigate the effects of the frustration on 
the spontaneous magnetization and other ground state 
properties, like spin stiffness. In Sec. IV we study the 
finite temperature properties: susceptibility, correlation 
length and specific heat. In the last Section we comment 
on the implications of our results. 

II. QHRT EQUATION FOR THE Ji - J2 MODEL 

We briefly summarize the basic steps followed for the 
derivation of the QHRT evolution equation. A detailed 
description of the method can be found in Ref»22i. Be- 
ing interested in a parameter region where antiferromag- 
netic correlations dominate, it is convenient to include 
in the Hamiltonian a staggered external magnetic field 
which directly couples to the Neel order parameter. The 
Hamiltonian is then written in terms of a reference part 
Hq, which in our case just coincides with the external 
field, and the interaction term Hj, chosen as the Hamil- 
tonian of the Ji — J2 model 

H = Ho + Hj = -h e's-^^R + Hj (2) 

R 

where g is the antiferromagnetic wavevector and Hj is 
defined in With our choice of Hq, the properties 

of the reference system are known. The method ap- 
plies to all interactions Hj, assumed bilinear in some 
operator p{r) ,i.e. Hj — ^ f dxdyp{x)w{x — y)p{y); in 
the present case p{r) represents local spin variables and 
the interaction in Fourier space is ^(k) = 2Ji[coskx + 
cos ky] + 4J2 cos kx cos ky. As a first step we write a for- 
mal perturbative expansion of the partition function of 
the model Z — Trexp(— /3iJ) in powers of the inter- 
action w^. Here we do not make use of any specific 
property of the reference system correlation functions 
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but instead the diagrammatic perturbative expansion is 
carried out in full generality. Remarkably, the result- 
ing perturbation series is formally equivalent to that of 
a classical system in (2 + 1) dimensions, with the peri- 
odic "time" variable now belonging to the interval (0, /3) 
and the interactions described by the classical two body 
potential wdrjt) = w{r)6{t)/ f3. Mapping on classical 
systems are frequently employed to describe quantum 
models at phase transitions, where long wavelength fluc- 
tuations always dominate. However, we stress here that 
our approach does not rely on any continuum limit and 
preserves the microscopic character of the Hamiltonian: 
actually, the formal expansion of the partition function 
leads to a quantum-to-classical correspondence which is 
exact, valid in the whole temperature range and not lim- 
ited to the critical regime. 

The thermodynamics of the corresponding classical 
system can be studied within the Hierarchical Reference 
Theory (HRT) formalism^^, successfully employed in the 
framework of classical statistical physics. HRT has been 
developed as a way to implement the basic ideas of the 



momentum space RG directly to microscopic classical 
Hamiltonians. The main goal of HRT is to describe 
how the inclusion of fluctuations over different length 
scales, affects the free energy of the model. For this pur- 
pose we define a sequence of auxiliary systems (named 
Q — systems) characterized by a cut-off dependent poten- 
tial w{k)Q which coincides with the physically relevant 
one w(k) for k > Q and vanishes elsewhere??. Introduc- 
ing a momentum cut-off in the interaction inhibits the 
antiferromagnetic fluctuations at long wavelengths, pre- 
empting the occurrence of phase transitions. When the 
cut-off Q is reduced, the properties of the system con- 
tinuously evolve between a mean field description, where 
fluctuations are absent, and those of the fully interacting 
system, where the fluctuations over all length scales are 
taken into account. Following this route, we obtain an 
exact set of differential equations expressing how the free 
energy a*^ and the correlation functions change when the 
cut-off Q is varied. In particular, the evolution of the 
Helmholtz free energy is given by 



ifk 



{ln(l-F<|^(k,c.„)iD(k)) 



2/3 ^ 

In [(1 + F<^-(k, u;„)?2;(k))(l -f F§-(k + g, uj^)w{\^ + g)) + i^^^(k, c^„)i^^^(k + g, c.„)^(k)^5(k + g)] } 



(3) 



r 



The summation is over the Matsubara frequencies LOn = 
27rn//3 and the integration is restricted to the one dimen- 
sional domain Sq defined by 



-I- C0S(fc.y)2 



(4) 



with Q G [0, 1]. Physically, the right hand side of this 
differential equation provides the contribution to the free 
energy due to fluctuations characterized by a wave vec- 
tor belonging to the domain Sq, which in fact spans the 
whole Brillouin zone when the cut-off parameter Q varies 
from 1 to 0. The choice of the sequence of surfaces is in 
principle arbitrary but definition proves quite conve- 
nient because it respects the symmetries of the lattice, 
guarantees the stability of the partial differential equa- 
tion © and leaves the critical fluctuations, of wave vector 
g, at the end of the integration. 

The key quantities appearing in Eq. ||2Jl are the spin- 
spin dynamical correlation functions of the Q-system 
i^^''(k,w„) =< 5"(k,w„)5'5(-k,-w„) >. When Q = 1 
no fluctuations have been included yet and the spin-spin 
correlation functions reduce to their expressions in mean- 
field approximation, obtained by keeping only the zero- 



loop contributions in the diagrammatic expansion: 

/i_L + w(k + g) 



i^Q=i(k,c.) 



m '^uj'^ + {fi± + w{'k + g)){^± + wCk)) 
m"2^2 ^ ^ ^(Ij. ^ g))(/i^ + w{\i)) 



Mil + w(k) 



(5) 



where /i_L and are known functions of temperature T 
and staggered magnetization m. For S = 1/2: /i^ = 
2rm-Hanh"^(2TO) and ^|| = 4T {1 - 4m'^)-\ The 
mean field form (O of the correlation functions corre- 
sponds, after Wick rotation, to the single-mode approxi- 
mation for the dynamic structure factor 



uj S{u — ei^) + S{lj + ek) 



(6) 



1 - e-^'^ /i_L + w{k) 
where the quasi-particle dispersion ek is given by 

Ck = m^/iwik) + fi±){w{k + g)+ (7) 

As pointed out in a better parameterization of the 
spin-spin dynamical correlation function Q is in princi- 
ple possible, allowing more complex dispersion relations. 
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Here, for simplicity and in order to verify the reliabil- 
ity of the method for frustrated systems, we consider 
the simpler approximation reproducing the lowest or- 
der in SWT. Linearization of Eq. (Q around k = 
shows that in mean field approximation, the spin veloc- 
ity Cs is proportional to the staggered magnetization m: 
Cs = m Ji -^8(1 — 2a), which correctly reproduces only 
the zeroth order SWT. 

For a generic value of the cut-off Q ^ 1, FQ^{k,uj) 
are not known because they already include the contri- 
bution of short range correlations. As a consequence, 
equation Q cannot be solved by straightforward integra- 
tion and additional information on the form of the dy- 
namical correlation functions is required. One possibility 
is to write an infinite set of exact differential equations 
also for the many particle correlation functions, thereby 
generating a full hierarchy of evolution equations: the 
QHRT. However, it is more convenient to limit our at- 
tention to the first equation of such a hierarchy, i.e. to 
Eq. ©, attempting a closure by imposing a suitable 
approximate parametrization of the spin-spin correlation 
function. The simplest choice, which we will adopt in 
this work, is to retain the same analytical structure of 
the mean-field approximation (jSJ keeping ^± and as 



free, cut-off dependent parameters. Their initial value at 
Q = 1 is set by the mean field expressions, while at Q ^ 1 
they are obtained by use of the exact sum rules which re- 
late the dynamical correlation functions to derivatives of 
the free energy a*^: 

^1 m = F^-(g,c. = 0) = [^2-4(Ji-J2)]-i (8) 



^] '=^^5^(g,^ = 0) = [/iff-4(Ji-J2)]-i (9) 

Within such an approximate closure, the inclusion of fluc- 
tuations does not change the analytical structure of Fq^: 
fluctuations lead to a renormalization of the parameters 
/ij and but do not allow for the creation of an inco- 
herent contribution to the dynamical structure factor 
By use of the representation jSJl in the exact evolution 
equation Q and summing over the Matsubara frequen- 
cies, we obtain an approximate, closed evolution equation 
for the free energy a*^ as a function of the cut-off Q and 
of the staggered magnetization m: 



dQ 2/3 



sinh 



sinh (^\f3m^{fi±_ + w(p))(^^ + w{p + g))^ 



-In 



+w'(p))(mii +u)(p-f-g)) 



(10) 



Equation H1Q(I . together with the thermodynamic consis- 
tency relations H8I9|) and the mean field initial condition 
at Q = 1, has been solved numerically for different val- 
ues of the temperature T and spin S. The Q-dependence 
of the parameters /i_L and ^|| is understood in Eq. (^UJ. 
We recall that Q is just a control parameter which al- 
lows for a gradual introduction of magnetic fiuctuations: 
All physical quantities we are interested in refer to the 
(5 = limit. However, the numerical determination of 
the free energy at Q = requires the integration of the 
parabolic partial differential equation in the whole inter- 
val < Q < 1. All the numerical results shown in the 
following sections will refer to Q = 0. In order to simplify 
the notation, the label Q will be dropped when we refer 
to Q = 0. We observe that the form of equation H10|) is 
unaffected by S and the spin value just enters through 
the initial condition at Q = 1. Equation H10|) has been 
written in quasi-linear form and solved by a fully im- 
plicit finite difference scheme. We stress that Eq. IjlUf) 
has been solved only at non-zero temperature and we did 
not attempt a direct analysis of the limiting equation as 
T ^ 0. 



III. ZERO TEMPERATURE LIMIT 

We first investigate which information on the ground 
state properties of the frustrated antiferromagnet may be 
extracted form the QHRT formalism. Let us consider the 
behavior of the zero field staggered magnetization when 
fluctuations are included. By differentiating Eq. H10|) with 
respect to the magnetization it is possible to obtain an 
evolution equation for the staggered field hq (m) , defined 
as the derivative of the free energy with respect to 
the staggered magnetization, at fixed m. This quantity 
implicitly describes the dependence of the staggered mag- 
netization mq on the cut-off Q at fixed magnetic field h 
through the identity hq{mq) — h. Note that in mean 
field approximation the model displays a finite tempera- 
ture critical point between a high temperature paramag- 



a = J2/J1 





0.1 


0.2 


0.3 


0.35 


0.4 


mo (S = 1/2) 


0.35 


0.34 


0.30 


0.25 


0.18 


0.05 


mo {S = 1) 


0.84 


0.81 


0.75 


0.64 


0.52 


0.34 



TABLE I: Ground state staggered magnetization mo for S = 
1/2 (top panel) and S = 1 (bottom panel). 
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FIG. 1: Q-evolution of the zero field staggered magnetization 
for S — 1/2 and different values of the frustration, at T = 
O.lJi. 
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■ 2 ■ S=l QHRT 
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0.1 0.2 0.3 0.4 0.5 

FIG. 2: Extrapolated ground state staggered magnetization 
mo divided by its classical value S. QHRT results for S = 1/2 
(empty circles) and 5 = 1 (full squares). First-order (empty 
triangles) and second-order (stars) SWT results from—. Lines 
are a guide for the eyes. 



net and an antiferromagnetic ground state. As a conse- 
quence, the initial condition of our QHRT equation shows 
non vanishing spontaneous magnetization toq=i ^ for 
/i = 0. According to the Mermin- Wagner theorem, ther- 
mal fluctuations depress, and eventually kill, long range 
Neel order and, as we know, QHRT does conform to this 
scenaric^°. Therefore, we expect that, at low temper- 
ature, the Q evolution of the staggered magnetization 
in zero field is characterized by a smooth initial behav- 
ior which extrapolates to the T = order parameter, 
followed by a sudden drop close to Q = due to the 
action of long wavelength fluctuations. This expected 
qualitative behavior can be checked on the numerical so- 
lution shown in Fig. Note that the QHRT equations 
preserve the convexity of the free energy in the physical 
limit Q = 0: This implies that at all temperatures T ^ 0, 
the zero-field inverse susceptibility must be strictly 
positive at Q = and then it has to cross Xq^ = at a 
non vanishing cut-off Q* 0. As a consequence, also the 
spontaneous magnetization, present at mean field level 
even at T 7^ 0, vanishes at Q* / as shown in Fig. 
n The main effect of frustration on the evolution of the 
spontaneous magnetization is in fact to move this charac- 
teristic cut-off toward larger values due to the expected 
decrease of the correlation length. 

Our estimates of the ground state staggered magne- 
tization, based on Fig. ^ are reported in Table U for 
both S = 1/2 and S = 1. The same data are shown 
in Fig. 121 together with the first and second order spin 
wave theory results^i. As already discussed in Refi^, at 
a = and 5=1/2 the QHRT estimate is very close to 
the commonly accepted value too ~ 0.307 However, 
while first order SWT predicts a vanishing magnetization 
at ac ~ 0.38, our results indicate that the ground state 
magnetization extrapolates to zero only at significantly 
larger values of a ~ 0.43. 

At zero temperature, in a symmetry broken phase, the 



spin stiffness ps can be defined in terms of the small mo- 
mentum behavior of the transverse, static structure fac- 
tor: 

S^k) - (11) 
2cs 

Within QHRT, the spin velocity Cg has been previously 
evaluated from Eq. (0 while S^_{k) is easily obtained 
from the parametrization |SJ|. As a result we obtain the 
simple analytical relationship: ps = TOq( Ji — 2J2) which, 
by use of Tab. Q provides an estimate of the stiffness. 
The proportionality between p^ and mg guarantees that 
at the quantum critical point, where the order parameter 
Too vanishes, the spin stiffness necessarily tends to zero, 
as expected. 

IV. TEMPERATURE DEPENDENCE OF 
PHYSICAL PROPERTIES 

According to the definition of the spin-spin correlation 
function, we can define the staggered magnetic suscepti- 
bility in terms of the zero frequency value of i^"'^(k,w) 
at k = g 

Xl = F-^(g,0)-(/i^-4(Ji- J2))-i 

Xfi = F-(g,0) = (M|| -4(Ji- J2))-i. (12) 

The results at to = are plotted in Fig. O as a func- 
tion of temperature for several values of the frustration 
ratio. In agreement with Mermin- Wagner theorem, true 
long-range order is absent at any finite temperature in 
two dimensions and in fact the susceptibility obtained 
via QHRT remains finite at every T ^ 0. Conversely, the 
divergence of x when zero temperature is approached is 
consistent with the onset of long range antiferromagnetic 
order in the ground state of the model. In order to better 
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the 5 = 1/2 system and different values of tfie frustration ra- 



tio. Straiglit lines are linear fits of the low temperature data, 
FIG. 3: S = 1/2 low temperature staggered susceptibility ^j^j^j^ ^j^^^jy ^jj^^ ^1^^ exponential law Qualitatively 
versus temperature, for various values of the frustration at ^-^-^^^ ^^^^^^^ ^^^^ ^^^^ obtained for S = 1. 
m — (xl = xf|)- Lines are guides for the eyes. 



analyze such a divergence we will refer to a closely related 
quantity: the antiferromagnetic correlation length, ^ de- 
fined as |kp|^^, where k = g + ikp is the location of 
the poles of the function F"'^(k, ~ 0). The same ^ es- 
sentially governs the exponential decay of the equal-time 
correlation function at T 7^ 0. Close to the critical value 
of the wave vector k = g, the two-point correlation func- 
tion at k = g + q and oj = takes the following form 



F"^(g-f-q,0) 



1 



/i-4(Ji- J2)+g2(Ji_2J2) 



(13) 



leading to a correlation length algebraically related to the 
susceptibility (|12|l and simply expressed as a function of 
the parameter fi: 



2 ^ ^^ - 4( Ji - J2) 
Ji - 2 J2 



riJi-2J2) 



(14) 



The numerical values of ^(T) can be easily deduced from 
Fig. 131 and clearly show a strong suppression of long 
range correlations when the frustration is increased. As 
predicted by weak-coupling (one-loop) RG analysis of the 
usual NLSM^^, when the temperature is lowered the sys- 
tem enters the classical renormalized regime character- 
ized by an exponential divergence of ^. In this region of 
the phase diagram the correlation length diverges expo- 
nentially, just as it would in the classical model: 



e ~ exp(2^p,/r) 



(15) 



where the coefficient ps is identified through the weak 
coupling RG analysis, with the spin-stiffness renormal- 
ized by quantum fluctuation3iS42L2&. As shown in Fig. 
the QHRT equation well reproduces, for different 
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• p from Eq. (15) 




■ SWT 




^ m'(Jj - 2 J^) 




III, ----^ 



FIG. 5: 5=1/2 renormalized spin stiffness as a function 
of the frustration. We report the values obtained from Eq. 
lllll (triangles), the coefficient pa from the low temperature 
divergence of the correlation length, Eq. II15II (circles) and 
the SWT estimate (squares) from—-. 



values of the frustration, the exponential behavior of ^ 
at low temperature as already found in Ref.'^". Eq. H15|) 
is valid in the low temperature limit (T << J) when cor- 
relations at the antiferromagnetic wave vector become 
long ranged. In this regime it is possible to estimate the 
renormalized stiffness coefficient which governs the ex- 
ponential divergence of the correlation length (|15() via a 
linear fit of the data in Fig. 0] 

The results of such a fit are shown in Fig. [3 together 
with the direct evaluation of the stiffness ps via Eq. 111|) 
at iS = 1/2. Frustration depresses the correlation length, 
nevertheless ps does not vanish in the range of frustration 
[0, 0.4] we have investigated and no clear sign of quantum 
critical point can be inferred from these data. A simi- 



lar result was actually found in Refr-- using variational 
spin-wave approach while exact diagonalization studies, 
followed by finite size scalingf2iS£ are not conclusive, the 
spin stiffness estimate depending on the chosen cluster 
geometry. 

The comparison between ps and the direct evaluation 
of the stiffness ps via Eq. (see Fig. O clearly shows 
that, within QHRT, the exponential divergence of the 
correlation length at low temperature is not governed by 
the same quantity which appears in the zero temperature 
transverse structure factor, contrary to the prediction of 
the RG analysis of the non linear sigma model^. In order 
to investigate the origin of such a discrepancy, we per- 
formed a careful analysis of the QHRT equation in the 
classical limit at J2 = 0, which corresponds to the S" — > 00 
case. At zero temperature, the ground state magnetiza- 
tion is clearly mo — S, providing the correct limiting 
result for the classical spin stiffness: ps = S^Ji. Instead, 
the numerical integration of the evolution equation gives 
a diverging correlation length consistent with Eq. H15(l 
where the renormalized stiffness is given by = 
This estimate actually coincides with the lowest order 
term of a large N expansion in the SU(iV) non linear 
sigma mode& evaluated at = 3. The correct RG re- 
sult for the SU(iV) non linear sigma model, valid at all 
orders in N, is instead given by ps = S'^Ji/{N — 2) which, 
at = 3 coincides with the zero temperature evaluation 
of Ps. A careful analysis of the RG approach in 2 + e 
dimension"^*^ suggests that the origin of the discrepancy 
found in the QHRT equations may be traced back to the 
chosen parametrization of the correlation functions 
The reason is twofold: 

• In the QHRT equation Hl()(l both longitudinal and 
transverse correlations contribute to one loop order, 
while in the RG approach only the transverse propagator 
enters the beta-function. Due to the non linear sigma 
model constraint, at low temperature the longitudinal 
correlations can be expressed in terms of the transverse 
ones: 



0.2 



< [l-7r(r)2]^[l-7r(0)2]3 > 
< [l-7r(0)2]5 >2 
i [< 7r(r)2^(0)2 > - < ^(0)2 >2] 



(16) 



where 7r(r) identifies transverse spin fluctuations. 
• Moreover, within our approximation Q the wave- 
vector dependence of the correlation function is not 
renormalized by fluctuations, while the RG analysis in 
2 + e dimensions shows that the one loop corrections do 
affect the momentum dependence of the one particle ir- 
reducible two point functior 
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FIG. 6: Plot of the function A(T) defined as A(T) = 
|log(C(r)) - lfit{T)\, where IfitiT) is the linear fit function 
describing the asymptotic behavior of log^(r) at low T. The 
function is plotted for different value of the frustration ratio. 



where A is the ultraviolet cut-off of the 0{N) non linear 
sigma model. This expression should be compared with 
the long wavelength limit of the corresponding QHRT 
approximation: 



r(p) = ppi_ + f3w{p + g) 



(18) 



showing that in QHRT the momentum dependence is 
fixed by the spin-spin interaction. 

These remarks suggest that a more flexible 
parametrization of the spin correlations should be 
used in the QHRT equations in order to achieve a 
quantitative representation of the low temperature 
divergence of the correlation length. 

In order to better understand the effects of frustration 
on the onset of Neel order, we also investigated the a de- 
pendence of the temperature Tx below which Eq. (|15|) is 
valid. In the following we refer to Tx as crossover temper- 
ature. At mean-field level we known that the inclusion 
of a n.n.n. interaction increases the bare value of the 
coupling constant g and tends to destroy the Neel order. 
When the system approaches the quantum critical point, 
controlled by the g = gc fixed point, the RG picture 
(qualitatively) predicts a decrease of the crossover tem- 
perature. The crossover temperature can be estimated 
by studying the temperature dependence of the function 
A(T) defined as A(T) = |log(C(r)) - lf^t{T)\, where 
lfit{T) is the linear fit function describing the asymp- 
totic behavior of log ^(T) at low T: A(r) then measures 
the deviations of £,{T) from the exponential law H15|l . 
Figure clearly shows that it is possible to identify a 
crossover temperature Tx by the requirement A(T) w 
for T < Tx. The results reported in Fig. Elseem to in- 
dicate that Tx does not depend much on a, its value 
remaining basically unchanged when frustration is in- 
creased. The weak dependence of the crossover tempera- 
ture on frustration may suggest either a very narrow crit- 
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FIG. 7: Zero field specific heat for S = 1/2 



ical region near the quantum transition, or a shift in the 
quantum critical point toward larger values of a > 0.5^^. 

Finally, from the free energy we calculated zero field 
specific heat C{T) for 5=1/2 and different frustrations, 
shown in Fig. Q. As expected, frustration play a mi- 
nor role at high temperature, the specific heat showing a 
paramagnetic behavior as 1/T^ for all values of the cou- 
pling J2. The position of the maximum moves at low 
temperatures and its intensity decreases for increasing 
value of the frustration, in agreement with SWT ^ and 
small-size Exact Diagonalization results—. 

V. CONCLUSIONS 

We have studied finite temperature properties of the 
frustrated Heisenberg model on a square lattice in the 
Neel phase, by use of the quantum hierarchical reference 
theory, a non-perturbative implementation of the renor- 
malization group method. We derived an exact evolu- 
tion equation which describes how the free energy of the 
system changes when antiferromagnetic fluctuations over 
growing wavelengths are included. This equation is sup- 
plemented by an approximate representation of the dy- 
namical correlation functions of the model Within 
the limitations induced by our parameterization of the 
two-point correlations, the method employed does not 
suffer from uncertainties in the extrapolation to the ther- 
modynamic limit^^. Our results are valid in a wide range 
of temperature and for all spin values because our ap- 
proach does not rely on any mapping on effective long 
wavelength theory and preserves the microscopic features 
of the model Hamiltonian. We focused on the role of fluc- 
tuations of wavevector close to g = (tt, tt) and then we 
could not investigate the possible occurrence of other or- 
der parameters characterized by different periodicities, as 
in spiral or coUinear phases. 

Starting from the knowledge of the free energy of the 



system, we calculated spin susceptibility, two point cor- 
relation function and correlation length at finite temper- 
ature and we also estimated ground state properties such 
as spin stiffness and staggered magnetization by extrap- 
olation to the T ^ limit. Calculations have been per- 
formed for spin S = 1/2 and S — 1 and different values 
of the frustration ratio, in order to investigate the role 
and extent of quantum corrections on top of the semi- 
classical spin-wave picture. The investigated range of 
frustration is restricted to a < 0.4 because beyond this 
point, fluctuations at wavevector (tt, 0) are expected to 
become important and our simple parametrization of the 
dynamical correlation functions cannot be retained. 

The spontaneous staggered magnetization extrapolates 
to zero at about a ~ 0.43 for S = 1/2 while a stronger 
frustration is required for S = 1. Within the QHRT ap- 
proach, both the spin velocity and stiffness are related 
to the order parameter and therefore are expected to 
vanish at the quantum critical point. For comparison 
we recall that in the non linear sigma model representa- 
tion of the quantum antiferromagnets, when the T = 
critical point is approached, both the staggered magne- 
tization and the spin stiffness vanish, while the spin ve- 
locity is expected to remain finite. Our results indicate 
that frustration strongly suppresses the correlations at fi- 
nite temperature but, at the same time, the value of the 
crossover temperature Tx remains basically unchanged 
when the frustration is increased. This fact may suggest 
a very narrow critical region near the quantum transition 
or a shift of the quantum critical point toward larger val- 
ues of a > 0.5. The estimate of the renormalized stiff- 
ness ps obtained via the low temperature divergence of 
the correlation length, Eq. (|15|l . does not vanish in the 
whole frustration range we investigated and is remark- 
ably larger than the direct evaluation of the stiffness ps 
obtained through the zero temperature static structure 
factor (|ll|l . The equality of these two quantities is com- 
monly expected on the basis of the known behavior of 
the SU(A^) sigma-model in the iV — > 00 limit and the 
RG analysis in 2 -I- e dimension—. The difference in the 
two estimates of ps we found in the SU(2) case is most 
likely a consequence of the approximate closure (jSJ we 
have studied, as shown by a detailed analysis of the clas- 
sical limit of the QHRT equations. In order to go beyond 
the approximation examined in this work, a representa- 
tion of the low temperature spin correlations richer than 
those commonly adopted in the many body literature^^ 
is needed. 

Nowadays, real materials displaying the properties of 
the Ji — J2 model in the Neel phase have not been yet 
synthesized. Nevertheless, the recent experiments at high 
pressure performed on the Li2VOSi04^ seem to indi- 
cate that a wide range of the Ji — J2 phase diagram 
can be studied by applying high hydrostatic pressure, 
which modifies the frustration ratio J2/J1 by changing 
the bond lengths and angles. Although the relationship 
between pressure and frustration ratio is not obvious at 
all, we believe that the results of our investigation may 
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be useful for the interpretation of future neutron scat- 
tering or NMR-NQR experiments. In addition ground 
state configuration of the VOM0O4 compound is still a 
controversial issue, first principle calculation of the J2 
coupling being in disagreement with the estimates based 
on the NMR-NQR experiments. Specific heat data 3* for 
the VOM0O4 compound may be analyzed on the basis of 



our finite-temperature results, the position of the maxi- 
mum in the zero-field specific heat giving an alternative 
way for estimating the frustration parameter. 

We thank L. Reatto and P. Carretta for useful discus- 
sions. The financial support of MIUR through PRIN- 
2003 is also acknowledged. 
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